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ABSTRACT 

The vertical warp in the debris disk /3 Pictoris - an inchned inner disk extending into a flat outer disk 
- has long been interpreted as the signpost of a planet on an inclined orbit. Direct images spanning 
2004-2010 have revealed /3 Pictoris b, a planet with a mass and orbital distance consistent with this 
picture. However, it was recently reported that the orbit of planet b is aligned with the flat outer 
disk, not the inclined inner disk, and thus lacks the inclination to warp the disk. We explore three 
scenarios for reconciling the apparent misalignment of the directly imaged planet /3 Pictoris b with 
the warped inner disk of {3 Pictoris: observational uncertainty, an additional planet, and damping of 
planet b's inclination. We find that, at the extremes of the uncertainties, the orbit of /3 Pictoris b has 
the inclination necessary to produce the observed warp. We also find that if planet b were aligned with 
the flat outer disk, it would prevent another planet from creating a warp with the observed properties; 
therefore planet b itself must be responsible for the warp. Finally, planet b's inclination could have 
been damped by dynamical friction and still produce the observed disk morphology, but the feasibility 
of damping depends on disk properties and the presence of other planets. More precise observations 
of the orbit of planet b and the position angle of the outer disk will allow us to distinguish between 
the first and third scenario. 

Subject headings: planet-disk interactions — stars: individual {j3 Pictoris) — planets and satellites: in- 
dividual (/3 Pictoris b) 



1. INTRODUCTION 

The (3 Pictoris debris disk is a rich system, with 
observational features resulting from the interplay of 
gravity, radiation pressure, collisions, infalling comets, 
sculpting by planets, and the physical properties of the 
gas, dust, and rocks that comprise t he disk. In the 
quart er-century following its discovery ([Smith fc Terril^ 
Il984f ). j3 Pictoris has epitomized young planetary sys- 
tems, amenable to state-of-the-art observations and to 
modeling of planetary formation processes. For a review, 
see lVidal-Madiar et al. (1998). 

A striking vertical warp in the /3 Pictoris disk, at 
approximately 85 AU from the sta r appears in op- 
tical and near-in f rared im ages (e.g. iHeap et al.l EoOffl 
Golimowski et al.' '20061). iMouillet et al.l ()1997l ) and 
Augcrcau et al. (2001) demonstrated that a perturbing 
planet on an inclined orbit could produce the warp, 
whose distance constrains the posited plan e t's mass and 
position. A decade later, iLagrange et al.l (|2009l I2010H 
discovered, via direct imaging, (3 Pictoris b, a planet 
consistent with producing the warp, if correctly inclined 
relative to t he disk. 

Ho wever. iCurrie et al.l (120111) - stitching to gether re- 
cent (jOuanz et al.ll2010l : iBonnefov et al.ll201lD data and 
newly reduced data (collected bv ILagrange et aril2010( ) 
- recently measured the planet's astrometric orbit and 
reported it to be, surprisingly, misaligned with the warp. 

lrdawson®cfa. harvard . edul 



They urged revisiting whether planet b could produce the 
warp and suggested an undiscovered additional planet as 
an alternative culprit. 

Here we explore three scenarios to reconcile the appar- 
ent misalignment of planet b with j3 Pictoris's warped 
inner disk. In Section [21 we model an inclined planet 
sculpting a planetesimal disk. Then (Section [3]), we con- 
sider the first scenario: within the extremes of the ob- 
servational uncertainties, /? Pictoris b has the inclination 
to produce the warp. In Section 21 we evaluate the pos- 
sibility that another planet created the warp; however, 
we find that the presence of planet b on a flat orbit pre- 
vents another object from creating the observed warp. In 
Section [5l we explore the third possibility: planet b had 
a higher inclination in the past, created the warp, and 
then its inclination damped. Thus we suggest (Section 
[6]) that planet b produced the warp, whether or not its 
orbit is currently aligned with the warped inner disk. 

2. MODEL OF A DEBRIS DISK SCULPTED BY AN 
INCLINED PLANET 

The warping of the (3 Pictoris debris disk by a planet 
on an inclined orbit results from the secular evolution 
of the disk's component planetesimals. The planet sec- 
ularly forces the disk, causing the planetesimals to os- 
cillate about the planet's inclined plane. Planetesimals 
in an initially flat disk (inclination z = relative to the 
plane of the flat outer disk) reach a maximum i of twice 
the planet's inclination ip. The oscillation frequency de- 
creases with the planetesimal's semi-major axis: in a 
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Fig. 1. — Snapshot (8 Myr) from an N-body simulation and an analytical model of an inclined planet sculpting an initially flat planetesimal 
disk. Top: Planetesimal inclinations vs. semi- major axis. The red lines indicate the planet's inclination ip and 2ip. The orange, blue, and 
black line is composed of triangles, each marking the inclination of a corresponding planetesimal in the middle panel. The planetesimals 
are color-coded: orange triangles have completed an oscillation, blue are just reaching their maximum inclination, and black are still at 
low inclinations. We use our analytical model (eqn. [ij to calculate the dashed gray line. Middle: Projected positions of the planetesimals. 
Lower: Same as middle using the analytical model; it reproduces the N-body simulation well. Compare these two lower panels to images 
of the disk, such as those shown in; Heap et a l. (2000), Fig. 8 and Golimowski et al. (2006), Fig. 5. 



young system like /3 Pictoris, the inner planetesimal disk 
- whicli secularly evolves quickly — is centered on the 
planet's inclined plane, while the outer disk - which sec- 
ularly evolves slowly - is still near its initial low incli- 
nation. Thus secular evolution produces an inclined in- 
ner disk aligned with the planet's inclined plane, a flat 
outer disk, and a warped feature, extending to 2ip, at 
the distance, Owarp, where planetesimals are just reaching 
their maximum inclination. This distance constrains the 



disk's evolution time and the perturbing planet's mass 
and location. In Fig. [H we plot the final inclinations 
and projected positions from an N-body simulation of 
an initially flat planetesimal disk sculpted by a planet 
on an inclined orbit. Each of 6000 test-particle plan- 
etesimals begins with e = i = 0, a semi- major axis a 
between 20-200 AU, and a random periapse, longitude 
of ascending node, and mean anomaly. The planet has 
the observed mass (mp — 9mjupitcr) and semi-major 
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axis (flp — 9.5 AU) of /3 Pictoris b, and ip = 3.6°; 
the star has mass m^, = 1.75msun. We used the Mer- 
cury 6.2 (|Chainber£lll999[ ) hybrid integrator with a step 
size of 200 days over a timespan of 8 Myr. The sys- 
tem's age, Vlt\ Myr (jZuckerman et aLll2001l), minus the 
planet formation time, < 3 — 5 Myr (|Hernandez et al.1 
I2007t iCurrie et alll2009[ ). imply that warp production 
likely began 3-20 Myr ago. We choose the viewing orien- 
tation to make both the outer disk and the planet's orbit 
edge-on, as observed. 

2.1. Planetesimal inclination evolution 

Secular evolution times set the warp's position: at the 
pointy part, planetesimals, in the midst of their first cycle 
about the planet's inclined plane, are just reaching their 
maximum 2ip (blue triangles. Fig. [1}. The components 
p and q of the planetesimal's instantaneous inclination, 
i = \/p^ + q-^, evolve as: 

ifroo sin(/i -I- 7) + 

Pforccd 

Zfree COs(/t -f 7) + ^forced 

n rUp 
4 m« 



P 



(1) 



where 



(2) 

The function 6 is a standard Laplace coefficient 
(jMurrav fc Dermottl[2000t) . and n = (Gm,/a3)i/2 ig the 
planetesimal's mean motion. When a single inclined 
planet forces an initially cold disk, the forced plane, with 

inclination iforccd = V Priced + ^Lccd' the inclined 
plane of the perturbing planet, ip, relative to the flat 
outer disk. Thus a planetesimal's inclination, initially at 
1 = 0, oscillates as: 



i = 2ip\ sm{ft/2)\ 



(3) 



2.2. Constraining the sculpting planet's orbit 

The disk reaches its maximal vertical extent, z^arp, 
at flwarp. From eqn. [3l the warp's secular oscillation 
frequency / is related to the disk's evolution time r by 
/t = vr, constraining the mass and semi-major axis of 
the perturbing planet: 

, _ "-warp rUp (1) 
" / — , '-*warp'-*warp^3 /2 V^warp y V 7 



and the perturbing planet's inclination is separately con- 
strained by: 



tan(2ip) 



(5) 



Our eqn. [4l is equ ivalent to the warp condit ion in 
iMouillet et alT (jlDQTl) , Section 5 and .Augcrcau et al.l 



(|2001f) . Section 1.2, de r ived f rom tidal theory. Addi- 
tionally, IMouillet et al.l (|1997[ ) modeled the disk's evo- 
lution using hydrodynamic simulations. However, using 
the secular equations above, we can not only constrain 
the parameters of a perturbing planet but produce the 
time-evolving disk morphology (e.g. Fig. [l] panel 3), 
facilitating comparisons to observations without simula- 
tions. 

The warp-revealing visual observations measure the 
light scattered by dust, likely produced by recent col- 
lisions of the planetesimal parent bodies. These dust 
grains are subject to radiative pressure, which induces 
larger, eccentric orbits relative to their parent bodies. 
Thus radiation pressure effectively increases the distance 
of the warp by a factor of ~ 1/(1 — 2f3), where f3 is the 
ratio of the radiati on forces to gravitational forces (see 
iChia ng et al.l [200^. For example, if /? = 0.2, a warped 
disk of parent bodies extending only to 50 AU creates an 
observed warp at 85 AU. In Fig. O we plot curves of Op 
vs. rUp that produce a warp at 85 AU (eqn. [4]) for /3 = 
and /3 = 0.2. Dust-generating collisions affect the par- 
ent bodies' orbits by damping ifrcc- However, the largest 
parent bodies in the coUisional cascade, which recently 
experienced their first collision and have not had Zfroc 
significantly damped, set the maximum vertical extent 
of the warp. 
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Fig. 2. — Three curves (solid) of (ap,mp) that produce a warp 
at 85 AU for (top to bottom) r = 3, 9, 17 Myr. A planet with 
properties above the lowest curves significantly impacts the warp 
evolution over the system lifetime. Equivalent curves for observed 
dust grains with radiative forcing parameter /3 = 0.2 are plotted in 
dashed grey. Increasing /9 decreases the warp distance in the plan- 
etesimal parent population, effectively shortening the warp propa- 
gation time. 



3. PLANET B POSSIBLY ALIGNED WITH INNER DISK 

The framework in Section [5] allows comparison of the 
model to the observed warp morphology. In creating 
Fig. [TJ we used the observed nip and Op. We selected 
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ip = 3.6° , corresponding (eqn. [5]) to a maxi mum vertical 
extent of -Zwar p = 11 AU at a warp = 8 5 AU ()Heap et alj 
[2OOOI Fig. 8: (Golimowski et al.ll2006l Fig. 5). We plot 
our simulations from the perspective such that both the 
outer disk and the planet's orbit are perfectly edge-on, as 
observed, which constrains the orientation of the warp. 
We find that the distribution of planetesimals matches 
t he disk shape well. 

iCurrie et al.l ()2011[ ) reported that planet b is aligned 
with the flat outer disk and misaligned with the inclined 
inner disk. Indeed, the difference between two separately 
measured angles, (a) the intersection of the planets or- 
bit w ith the sky plane, 31.32°[30.56, 32.12] (|Currie et al.l 
I2011D. and (b) the sky position of the o uter disk, 30-31° 
(iKalas fc Jewit^ [1991 or 29.5° ± 0.5 f Boccaletti eTal] 
1200911 . is consistent with alignment. However, at ~ 2-cr, 
/3 Pictoris b may be misaligned with the flat outer disk 
by ip = 3.6°, and aligned with the middle of the inclined 
inner disk, as in our model (Fig. [T|). 

4. PLANET B PREVENTS ANOTHER PLANET FROM 
CREATING THE WARP 

If, contrary to the scenario of Section [Sj more precise 
measurements confirm that planet b's orbit is aligned 
with the flat outer disk, an undiscovered planet c would 
be an obvious suspect. However, the presence of planet b 
severely restricts the parameter space for an additional, 
warp-making planet. The following loose, generous con- 
straints on planet c must be simultaneously satisfied: 

1. To have escaped radial velocity detection by 
iGalland et al.l (|2006[ ). planet c must satisfy 

""^ < 9Vac/lAU 

^Jupiter 

2. Planets b and c must be sufficiently separated for 
stability, obeying 

— >2.4(M, + Mc)^/' 
a 

where fj, is the planet-to-star mass-ratio ([GladmarJ 
119931 ). Even if planets b and c underwent scatter- 
ing, it is unlikely that their unstable configuration 
would last long enough to create the warp. 

3. Planet c must create the observed disk morphol- 
ogy - an inclined inner disk from 40-90 AU - with- 
out exciting planet b's inclination. We generously 
allow any parameters for planet c that produce 
a forced inclination iforccd < 2° for planet b and 
3° < jforced < 8° for the inner disk. We calculate 
the for ced inclinations usiiig mu lti-planet secular 
theory (jMurrav fc Dermottll2000l . Section 7.4). 

4. Given the sy s tem age of 8-20 Myr 
(jZuckerman et all 1200 1[ ). planet c must have 
a mass and semi-major axis small enough so that 
the evolution time of the warp at 85 AU is slower 
than 1 Myr (but faster than 20 Myr). If the 
observed dust grains have /3 > (Section 12. 2p , the 
evolution time must satisfy this constraint at 85 
AU (1 — 2/3), a more restrictive lower limit. The 
lower limit is generous, allowing for the possibility 
that planet c became inclined very recently. 



5. A secular resonance should not disrupt the flat 
outer disk (90-200 AU); we calculate the locations 
of secular resonances using multi-planet secular 
theory. 

We explore the parameter space of planet c's semi- 
major axis, mass, and inclination. For each combination, 
we evaluate the equations for the five criteria above in 
a two-planet system containing planets c, and b, with 
its nominal nib = 9mjupiter, Ofc = 9.5 AU, and ib = 0. 
In Fig. [21 we shade the (oc, rric) regions for which no 
possible inclination of planet c can satisfy the constraints. 
Other choices for rub and af,, within observational errors, 
yield qualitatively similar results. Constraint 3 is most 
restrictive: the planet must have high enough mass to 
excite the warp, but low enough mass not to excite ib] 
it must be far enough out that the warp can extend to 
85 AU, but close enough to incline the inner disk at 40 
AU. Constraint 4 considers the time dependence: the 
warp must reach 85 AU in the system age. The secular 
oscillation frequency at 85 AU depends on the mass and 
semi-major axis of both planets b and c - but not on 
their inclinations. Planet b, even on a non-inclined orbit, 
makes a large contribution to this frequency, leaving little 
room for a contribution from planet c. 




Fig. 3. — Constraints on and rric- The region shaded 
horizontal-striped green violates Constraint 1 (lack of RV de- 
tection), upward-slanted black violates Constraint 2 (stability), 
downward-slanted red violates Constraint 3 (produces disk mor- 
phology without exciting planet b), vertical-striped blue violates 
Constraint 4 (timescale consistency), and shallow-slant purple vio- 
lates Constraint 5 (secular resonances in the outer disk). See text 
for details. 

Therefore we can rule out that an undiscovered planet c 
causes the warp, because it cannot do so in the presence 
of planet b. Additional planets may be present in the 
system, but are unlikely to be predominantly responsible 
for the warp. 

5. PLANET B'S INCLINATION MAY HAVE DAMPED 
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We demonstrated that planet b's current orbit is con- 
sistent with producing the warp only at the extremes of 
the uncertainties (Section[3l) and that planet b's presence 
inhibits another planet from creating the warp (Section 
3]). If follow-up observations confirm the nominal orbit of 
(5 Pictoris b and position angle of the fiat outer disk, we 
are left uncomfortably with a planet misaligned with the 
warp it produced. However, we need not regard that con- 
clusion with such discomfort. A sculpted planetesimal 
disk can record the history of a planet's orbit, revealing 
a dynamical past we would never guess from the planet's 
current orbit. For example, Neptune has a nearly circu- 
lar orbit today but may have have sculpted the Kuiper 
Belt, our solar system's remnant planetesimal disk, dur- 
ing a period of high eccentricity, whi ch was su bsequently 
damped by dynamical friction (see ILevison e t al. 2008, 
and references therein). Embedded in a planetesimal 
disk, P Pictoris b would experience damping of its incli- 
nation, though the extent and timescale depend on disk 
properties. 

5.1. Consistency with disk morphology 

First we demonstrate that the disk morphology can 
be consistent with the damping of /3 Pictoris b's orbital 
inclination. The planetesimals, with initial forced incli- 
nations of ip, begin to oscillate about the fo r ced p lane 
and, as we showed in lDawson fc Murrav-Clavl ()201lD . ef- 
fectively freeze at the inclination values they reach after 
one damping timescale of the planet's inclination. There- 
fore, the warp freezes at the distance it reaches when the 
planet's inclination damps. Because their forced incli- 
nations are damped, the planetesimals have a maximum 
i = ip instead of 2ip, requiring /3 Pictoris b to have an 
initial ip ^ 7°. Fig. SI left panel shows an example: 
the disk morphology still matches the observations even 
though the inclined inner disk is not aligned with the 
planet's orbit. 

However, if the disk evolves for too long after the 
planet's inclination damps, the planetesimal nodes ran- 
domize and the morphology becomes boxy (right panel of 
Fig. 01). The distinction between a consistent and an ex- 
cessively boxy morphology is qualitative and may require 
sophisticated modeling of the observed scattered light. 
However, it is clear that no more than a partial preces- 
sion period at 85 AU can have passed since the planet's 
inclination damped. Since the system is young, this re- 
quirement does not demand implausible fine-tuning. 

5.2. Damping conditions 

We consider the feasibility of damping planet b's or- 
bital inclination. We expect that planet b started in the 
outer disk's plane, was perturbed, and is returning to its 
original plane. Possible per turbations include a sec ond 
planet scattered inward (e.g. lJuric fc Tremain j|2008() , or 
reso nant-induced inclina tions in the disk near the planet 
(e.g. iThommes &: Lissau cr 2003). 

Ejection of a scattered planet could alter the total an- 
gular momentum, requiring a disk with mass comparable 
to planet b to keep the forced plane, to which planet b 
damps, low (e.g. an unusually long-lived gaseous proto- 
planetary disk - iHillenbrand et alTl 19931 - or a particu- 
larly massive planetesimal belt). We confirmed that fast 
precession of planet b caused by a massive disk, since 



the precession occurs about a misaligned axis, would not 
prevent the excitation of the warp. On longer timescales, 
the damping medium must allow planet b to remain in- 
clined for almost the system's age. High-mass planets 
whose inclinations bring them several scale-heights out 
of the disk fail to open a gap, and t hus gas interaction 
damp s them after only 100 orbits ()Marzari fc NelsonI 
l2009t) . A more moderate inclination (i < 10°) for planet 
b, allowing it to maintain a clea n gap, would lead t o 
satisfactorily- inefficient damping (jBitsch fc KlevI [20111) . 
perhaps on a 10 Myr timescale. 

However, a modest disk mass may suffice. During 
planet-planet scattering, the total angular momentum 
is conserved, producing an average forced-plane still 
aligned with the flat outer disk. Planetesimals contribut- 
ing to dynamical friction follow the forced plane, with 
each planetesimal attempting to damp the planet to the 
planetesimal's own plane. Even if planet b dominates 
the forced plane in its immediate vicinity, distant plan- 
etesimals, within a factor of several of the planet's semi- 
major axis, could damp the planet's orbit to their flat 
forced plane (i.e. the plane of the outer disk). (How- 
ever, planet b would need to dominate the forced plane 
from 4 - 90 AU to excite th e inclined inner disk.) Fol- 
lowing l^rT&IChianJ (|2003), Section 2.3, a surface den- 
sity as low as 3 lunar masses within 40 AU could damp 
ip from 10° over 4 million years if the disk remains 
thin. Collisional dis sipation could keep the disk thin 
(|Goldreich et al.|[2004l eqn. 33, 50) if 

s < idisk > mdisk ,„s 
< (6) 

1cm 2° 5mEarth 

where s is a typical planetesimal radius. Since the colli- 
sional dissipation rate is oc cr oc a^^, where a is the disk 
surface density, collisional dissipation could keep the disk 
thin near planet b while allowing excitation from 40-90 
AU, where the warp is created. 

Clearly further study is required to flnd evolutionary 
scenarios that produce the inclination damping used in 
Fig. m These details are independent of this section's 
main message: a transient planetary orbit could establish 
a warp in a disk, to be observed at present, even if that 
planetary orbit has since changed. 

6. CONCLUSION 

We explored three scenarios for the apparent misalign- 
ment between the warped inner disk of (3 Pictoris and 
the orbit of the directly imaged planet b. In the first, 
most plausible scenario (Section planet b's orbit is 
consistent with producing the warp, at the extremes of 
the uncertainties. We argued that the alignment depends 
not only on the planet's orbit but on the (separately mea- 
sured) position angle of the outer disk. Therefore both 
of these quantities must be measured more precisely and, 
if possible, simultaneously from the same images. 

In the second, most obvious scenario (Section 14]), an- 
other planet warps the disk. However, we demonstrated 
that planet b inhibits another planet from producing the 
warp. Other planets may exist in the system, creating 
other disk features, but they cannot be responsible for 
the warp. 

If the first scenario is ruled out by more precise ob- 
servations, we are left with the third scenario (Section 
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Fig. 4. — Snapshots at two times (left and right) from an N-body simulation of a planet (rrip = 12mjupiter; 
) = 6°) sculpting an initially flat planetesimal disk. We imposed damping of the planet's inclination of the form | = 



13.5 AU; initial 



27r/(4 Myr), following 

Appendix A of Wolff et al. I 120111) . Top: Planetesimal inclinations vs. semi-major axis. The red dashed and solid lines are the planet's 
initial and current inclination, respectively. Bottom: Projected positions of the planetesimals. See Fig. [T]for color coding. 



[5]): planet b created the warp and then had its inchna- 
tion damped. Detailed modeling of scenarios that allow 
for the damping of planet b's inclination will be neces- 
sary. Confirmation of the damping scenario, especially if 
observers discover more systems with planets misaligned 
with the warp they produced, could shed light on disk 
properties that are important for planet formation but 
difficult to measure directly. 
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